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Abstract 

North-east India, being a corridor to South-east Asia, is believed to play an important role in transmitting drug resistant 
Plasmodium falciparum malaria to India and South Asia. North-east India was the first place in India to record the emergence 
of drug resistance to chloroquine as well as sulphadoxine/pyrimethamine. Presently chloroquine resistance is widespread all 
over the North-east India and resistance to other anti-malarials is increasing. In this study both in vivo therapeutic efficacy 
and molecular assays were used to screen the spectrum of drug resistance to chloroquine and sulphadoxine/pyrimethamine 
in the circulating P. falciparum strains. A total of 220 P. falciparum positives subjects were enrolled in the study for 
therapeutic assessment of chloroquine and sulphadoxine/pyrimethamine and assessment of point mutations conferring 
resistances to these drugs were carried out by genotyping the isolates following standard methods. Overall clinical failures 
in sulphadoxine/pyrimethamine and chloroquine were found 12.6 and 69.5% respectively, while overall treatment failures 
recorded were 13.7 and 81.5% in the two arms. Nearly all (99.0%) the isolates had mutant pfcrt genotype (76T), while 68% 
had mutant pfmdrA genotype (86Y). Mutation in dhps AYJ codon was the most prevalent one while dhfr codon 108 showed 
100% mutation. A total of 23 unique haplotypes at the dhps locus and 7 at dhfr locus were found while dhps-dhfr combined 
loci revealed 49 unique haplotypes. Prevalence of double, triple and quadruple mutations were common while 1 haplotype 
was found with all five mutated codons ( F/AGEGS/T ) at dhps locus. Detection of quadruple mutants {51I/59R/108N/164L) in 
the present study, earlier recorded from Car Nicobar Island, India only, indicates the presence of high levels of resistance to 
sulphadoxine/pyrimethamine in north-east India. Associations between resistant haplotypes and the clinical outcomes and 
emerging resistance in sulphadoxine/pyrimethamine in relation to the efficacy of the currently used artemisinin 
combination therapy are discussed. 
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Introduction 

Control of Plasmodium falciparum (Pf) malaria in tlie north-east 
India, reporting 0.15-0.2 million cases of malaria and over 100 
deaths annually [1] is becoming challenging due to the widespread 
resistance to chloroquine (CQ) and increasing therapeutic failure 
to sulphadoxine/pyrimethamine (SP) [2]. Chloroquine resistance 
in Pf first appeared more or less simultaneously in Southeast Asia 
(Thai-Cambodian border) and South America (Colombia) in the 
late 1950s [3,4]. In India, resistance to chloroquine in Pf was first 
detected in 1973 in Diphu area of Karbi-Anglong district of Assam 



state, north-east India [5] which gradually spread towards the east, 
west and south, covering almost entire country [6]. While 
resistance to SP was first reported again from Thai-Cambodian 
border in 1960s [7], initial reports of SP resistance in India, 
emerged in 1979, again in Karbi-Anglong district, Assam [8] and 
later documented in 1992 in Changlang district of Arunachal 
Pradesh, north-east India [9] indicating epidemiological impor- 
tance of this region in spreading the drug resistant strains. 
Recently, treatment failure to CQ, SP and quinine, administered 
sequentially, was noted in 6% of Pf cases, indicating the presence 
of multidrug resistant Pf strains [10]. It has been observed that 
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though the slide positivity rates (SPR) and slide falciparum rates 
(SFR) in north-east India have increased following emergence of 
drug resistance, the trend of progression of SFR over time has 
remained more or less linear, the case fatality rate (CFR) showed a 
steady increase (Figure SI & Figure S2 as per north-east Indian 
states NVBDCP data). 

Resistance to anti-folates (e.g. sulphadoxine/pyrimethamine) 
has been found associated with the point mutations at dihydro- 
folate reductase (dhfr) and dihydropteroate synthetase (dhps) genes 
and that of chloroquine with Pf chloroquine resistance transporter 
{Pfcrt) and Pf multidrug resistance {Pfmdrl) genes, [11,12,13]. 
Mutations at several amino acid positions of the Pfcrt gene from 
the chloroquine resistant parasite lines have been reported, of 
which mutation at 76 po.sition from lysine to threonine; (K76T) was 
consistendy present among CQ;resistant P. falciparum isolates 
[14], making it a good marker to provide information on the Pf 
chloroquine resistance status. Additionally, point mutations at 
N86Y, Y184F, N1042D and D1246Y loci ot pfmdr-l gene are 
reported to reduce susceptibility to chloroquine [1,5,16]. Pyri- 
methamine inhibits the parasite enzyme Pfdhfr involved in folate 
biosynthesis pathway [17,18], However, parasites with point 
mutations at codon position 51, 59, 108 and 164 of dhfr gene 
require higher amount of pyrimethamine to inhibit its activity and 
subsequendy become resistant to pyrimethamine. The first 
mutation was found to be in position 108 (S108N). The drug 
sensitivity is reduced further if this mutation occurs along with 
mutations at N51I (Asn-Ile), C59R (Cys-Arg) or I164L (Ile-Leu). It 
was noted that the double mutation in dhfr gene (N51I-hS108N or 
C59R-I-S108N) causes intermediate resistance while the triple 
mutation of the codons S108N-I-N51I-I-C59R confers high 
resistance [19]. However, an additional mutation at codon 
I164L leading to a quadruple mutation results in highest level of 
resistance to SP as well as to chloroproguanil dapsone and 
artesunate dapsone proguanU [20]. Sulphadoxine competes with 
the substrate that binds to the parasite enzyme Pfdhps. Mutations 
at several amino acid positions of this enzyme have shown 
reduction in its binding capacity to the drug. Mutations in pfdhps 
gene have been reported at five different amino acid positions 
(436, 437, 540, 580 and 613). In dhps mutation first starts at codon 
S436F/A (Ser-Ala/Phe) or at A437G (Ala-Gly) followed by 
mutations at other positions viz., K540E (Lys-Glu), A581G (Ala- 
Gly) and A613S/T (Ala-Ser/Thr) conferring more resistance 
[1,21,22,23,24,25,26]. The presence of mutations at codons 437 
and 540 of pfdhps together with the triple mutation of pfdhfr 
[quintuple mutation] is a significant predictor of SP treatment 
failure [27], 

CQwas the drug of choice in north-east India, except in some 
select areas where SP was used as the first line of treatment till 
2008. Thorough assessment of the extent of anti-malarial drug 
resistance and therapeutic failure is required before taking a call 
on making any changes in the drug policy, particularly in highly 
endemic areas such as north-east India. In view of this, the present 
study was conducted during 2006-2007 in two districts of Assam 
state in north-east India covering two sites located at considerable 
distances to determine the therapeutic efficacy of CQ, and SP 
against uncomplicated falciparum malaria. Further, the frequency 
of mutations associated with CQ, and SP resistance were 
determined in circulating Pf strains to correlate their association 
with therapeutic efiicacy. 



Materials and Methods 

Study sites and study design 

The study was carried out in two malaria endemic districts of 
north-east Indian state of Assam viz., North Lakhimpur and 
Hailakandi (Figure 1). Fever cases attending the Primary Health 
Centre of these two study sites were screened for malaria infection 
using conventional microscopy by examining both thick and thin 
blood smears stained with 3% Giemsa. Uncomplicated falciparum 
malaria cases meeting the following inclusion criteria were 
included in the study (1) >6 months of age (2) initial parasite 
count within 1000-1,00,000 parasites/|xl on Day 0 (3) not suffering 
from severe malnutrition (4) mono infection with P. falciparum (5) 
no signs of any other illness and other causes of febrile disease (6) 
axillary temperature S37.5°C (7) written informed consent for 
willingness to participate in th(; study (8) no history of allergic 
reactions to sulfa-containing drugs and (9) able to return to clinic 
for follow-up visits. 

A total of 220 Pf positives subjects fulfilling above criteria were 
enrolled in the study at PHC Kathlicherra in Hailakandi district 
(n= 120) and PHC Nowboicha in North Lakhimpur district 
(n= 100) and the therapeutic assessment of CQ and SP was 
carried out. Both drugs were supplied by the Directorate of 
National Vector Borne Diseases Control Programme, Govern- 
ment of India. CQ was administered orally at 2,5 mg/kg body 
weight in three divided doses (10 mg/kg on Day 0, 10 mg/kg on 
Day 1 and 5 mg/kg on Day 2). SP was given orally as a single dose 
at 25 mg sulphadoxine and 1.25 mg pyrimethamine/kg body 
weight. Antimalarials were administered under direct supervision 
of the investigators and patients were foUowed-up for a period of 
28 days for CQand 42 days for SP. Patients showing failure were 
treated with parental quinine in 5% dextrose until they were able 
to take the drug orally. All patients were given 0.75 mg 
primaquine/kg body weight as a single dose, for radical cure, at 
the end of the study. 

Clinical and Parasitological evaluation 

The standard WHO protocol [28] for classification of response 
to anti-malarial treatment in low to moderate transmission areas 
was adopted. Blood smears were taken for parasite density count 
on day 0, 1, 2, 3, 7, 14, 21, 28 and 42 of recruitment along with 
the recording of axillary temperature and chnical assessment of the 
enrolled subjects. In addition, blood samples were collected from 
the willing subjects on filter paper on day 0 for molecular studies. 
The clinical response was classified as (i) early treatment failure 
(ETF) (ii) late clinical failure (LCF) (iii) late parasitological failure 
(LPF) and (iv) adequate clinical and parasitological response 
(ACPR). 

Sample processing and DNA Isolation 

Approximately 50 jo^l of blood from each fever case at the time 
of initial recruitment was spotted on Whatman no. 1 filter paper, 
desiccated and transferred to the laboratory'. Parasite DNA was 
extracted from 4-5 punches (3 mm diameter) cut from dried blood 
spot from each filter paper using QIAamp DNA Blood Mini Kit 
(Cat. 51106, QIAGEN, USA) following manufacturer's instruc- 
tions. The parasite DNA was eluted in 50 |a.l volume and stored at 
4°C for immediate use. The parasite species was determined using 
primers and method described earlier [29]. 

Parasite genotyping 

P. falciparum mutations associated with CQ and SP resistance 
were typed by nested PCR-RFLP method as described and 
available on the web (http://medschool.umaryland.edu/CVD/ 
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Hailakandi 



Figure 1. Map of the studied areas of north east India. 
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plowe.html). In this study, the following codons and polymor- 
phisms were analyzed: pfmdr 86, pfcrt 76; pfdhfr (16, 50, 51, 59, 
108, 164) and pfdhps (436, 437, 540, 581, 613). All PGR 
amplifications were performed with MJ Mini thermal cycler (Bio- 
Red) in a reaction volume of 25 |J.l with Promega master mix. 5 |a.l 
of extracted DNA was used as a template in 1 st round PGR while 
1-2 10,1 of primary PGR product was used as a template for nested 
PGR. 5-8 |j.l of 2nd round PGR product was digested with 1 unit 
of respective restriction enzyme (NEB) in a 20 |J,1 reaction volume 
using NEB buffer. Amplicons and fragments were separated on 
2% or 3% agarose gels, stained with ethidium bromide and 
visualized under UV. The detailed description of sample size used 
for different aspects are outlined in Figure 2. 

Data analyses 

The data was double entered in SPSS VI 3.0. All statistical 
analyses were carried out using SPSS V13.0 and XLSTAT V- 
2012.5.02. 

Mutations were determined on the basis of the presence or 
absence of the expected bands on ethidium bromide stained 
agarose gel analysis from PGR amplicon after restriction 
digestions. Site specific distribution of allele was expressed in 
percentage, dhps and dhfr haplotypes were counted separately and 
in combination and their prevalence recorded. Associations 
between dhfr/dlips alleles and SP therapeutic failure and pfcrt 
alleles and GQ, therapeutic failure were assessed using an adjusted 
Odds ratios (OR) and the 95% confidence intervals. Appropriate 
statistical tests were carried out. 

Ethical consideration 

The study was approved by the Scientific Advisory Gommittee 
and Human Ethical Gommittee of Regional Medical Research 
Gentre, Dibrugarh, Assam (India). All participants provided 
informed written consent. However, in case of a minor or children 
enrolled in the study, informed written consent was taken either 
from their parents or from their legal guardians, if parents were 
not available. As proposed by the Institutional Ethics Gommittee, 
Regional Medical Research Gentre, Dibrugarh, Assam, all 



participants were informed about the objectives and intended 
use of the findings prior to their enrollment in the study as well as 
their right to withdraw from the study at any time without 
notification. The informed consent was translated into local 
language. 

Results 

Therapeutic efficacy of CQ and SP 

A total of 200, among the 220 enrolled subjects with acute 
uncomplicated malaria, from the two sites were included in the 
study (50 for GQ_and 50 for SP from each site) and treated with SP 
and CQ. Patients' characteristics, demographic data and treat- 
ment outcomes are shown in Table 1 . Of these, 95 subjects each 
successfully completed the 28 days foUow up for CQ_ and 42 days 
follow up for SP. The therapeutic efficacy study at the two sites 
revealed 86.3% and 19% ACPR response with SP and GQ 
respectively. Overall clinical failures in SP and GQarms were 12.6 
and 69.5% respectively, while overall treatment failures were 13.7 
and 81.5% in the two arms. Hailakandi recorded the higher GQ 
failure while both the sites had almost similar rate of SP failure 
(Table 1). 

Prevalence of pfcrt and pfmdr^ mutations 

The status oi pfcrl K76T and pfmdr-\ N86Y mutations was 
determined in 200 (119 from Hailakandi, 81 from North 
Lakhimpur) and 196 (Hailakandi = 1 19, North Lakhimpur = 77) 
isolates from the two locations. Based on pfcrl locus almost all 
isolates were found resistant to GQ at both the sites [North 
Lakhimpur 80/81 (98.8%), Hailakandi 118/119 (99.2%)], while at 
the pfmdr-l locus 80% (61/77) isolates from North Lakimpur and 
61% (72/119) isolates from Hailakandi were found resistant. 
Overall, 99% isolates had mutant pfcrt genotype (76T), while 68% 
had mutant /)/miir- 1 genotype (86Y) locus (Figure 3). 

Prevalence of dhps and dhfr mutations 

The presence of five dhps mutations (S436F/A, A437G, K540E, 
A581G, and A613S/T) and six dhfr mutations (A16V, G50R 
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Genotype at pfdhps locus 




HLY: District IHailalcandi, Assam 
NLK: District North Laldiimpur, Assam 



Total 173 samples 



Figure 2. Sample size, site and genotyped codon wise. 

doi:l 0.1 371/joumal.pone.0105562.g002 



N51I, C59R, S1()8N and I164L) known to be involved in 
resistance to SP were evaluated (Table 2). Mutation analysis was 
successful in 173 subjects for dhps codons and in 168 patients for 
dhfr codons. 

Mutant codons were quite common at both the loci and study 
sites. Among all mutations studied, mutation in dhps 437 codon 
was the most prevalent, while mutation in dhps 581 codon was the 
least prevalent. Mutant alleles of dhps were more common in 
Hailakandi. The mutant dhfr 50 allele was completely absent 
while dhfr codon 108 showed 100% mutation (Table 2, Figure 4). 

Distributions of dhps and dhfr haplotypes 

The allelic pattern or haplotypes were constructed for dhps and 
dhfr loci individually as well as in combination. A total of 23 
unique haplotypes at the dhps locus (Table 3) and 7 haplotypes at 
dhfr locus (Table 4) were found. In case oidhps, 40 isolates (23.1 "b 
of isolates) had all five sensitive codons (SAKAA). On the other 
hand 5 haplotypes (22.5'!^)) with single mutation (SAKA S/T , 
SAKGA, SAEAA, SGKAA, F/AAKAA), 7 haplotypes (25.4%) 
with double mutations (SAK GS/T , SAEAS/T, SGKAS/T, 
SGKGA, SGEAA, F/AAEAA, F/AG KAA), 6 haplotypes 
(22.5%) widi triple mutations (SGK GS/T , SGEAS/T, S GEGA , 
F/AAEAS/T, F/AG KAA, F/AGEA A), 3 haplotypes (5.8'X,) 
with quadruple mutations (S GEGS/T , F/AGEAS/T , F/ 
AGEGA) and 1 haplotype (0.6%) with all five mutated codons 
(F /AGEGS/T) were found. Higher prevalence of resistant allele 



at codon 437 (14 haplotypes), codon 540 (12 haplotypes) and 
codon 613 (11 haplotypes) was noted suggesting that mutations 
first started at 437 codon position of dhps locus in NE India and 
subsequenfly associated with other mutations to confer resistance 
(Table 3). 

At the dhfr locus, only 7 haplotypes were observed among 168 
isolate (Table 4). No haplotype having all the 6 sensitive codons 
was observed. There was only 1 haplotype (1.8% of isolates) 
carrying single mutation (ACNCNI), 2 haplotypes (65.5%) with 
double mutations (ACNRNI, ACNCNL), 3 haplotypes (14.3%) 
with triple mutations (VCNRNI, ACIRNI and ACN RNL) and 1 
haplotype (1.8%) with quadruple mutation (AC IRNL) . The 
prevalence of haplotype ACN RN I containing mutations at codons 
59 and 108 had the highest frequency (64.3%). 

Construction ot dhps-dhfr combined haplotypes among the 154 
isolates, where PCR-RFLP results for all codons were available, 
revealed 49 unique haplotypes, of which SAKAA-ACNRNI 
(13.6%), SGKAA-ACNRNI (8.4%) and F/AGEA.4-ACNRNI 
(7.1%) haplotypes had the higher frequencies (Table 5). Assess- 
ment of dhps /dhfr combined haplotypes with single or multiple 
mutations viz., S108N, N51I+S108N, C59R+S108N, N51I+ 
C59R+S108N, C59R+S108N+I164L, N51I+C59R+S108N+ 
I164L, N51I+C59R+S108N/A437G, N51I+C59R+S108N/ 
K540E and N51I+C59R+S108N/A437G+K540E, considered as 
predictive for SP resistance or failure [26], revealed that C59R+ 
S108N (n = 108; Hailakandi 65%, Nortii Lakhimpur 59%), N51I+ 
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drug resistance molecular marker 

Figure 3. Prevalence of point mutations at pfcrt and pfmdrl 
locus in the two studied sites. 

doi:10.1371/journal.pone.0105562.g003 

C59R+S108N (n = 31; Hailakandi 1 7.5%, North Lakhimpur 19%) 
and C59R+S108N+I164L (n = 19; Hailakandi 6%, North 
Lakhimpur 19%) combinations were the most prevalent ones 
(Figure 5). 



Association between dhfr-dhps point mutations and 
treatment failure with SP and CQ 

SP treatment failure was strongly associated with dhps-dhfr 
mutations and their combinations. It was found that 5 dhps-dhfr 
combined haplotypes were associated with 6 cases of early 
treatment failures of SP, while 25 dhps-dhfr combined haplotypes 
were associated with 33 cases of early treatment failure of CQ,. 
Most of the treatment failure haplotypes were with mutations 
59R+108N or 51I+59R+108N or 59R+108N+164L at the dhfr 
locus indicating their role in SP treatment failure in north-east 
India (Table 6). Among the CQ failure cases, 67% belonged to 
adult age group (>15 years), 24% to 5-15 years age and 9% were 
under 5 years of age. However, in case of SP treatment failures, 
50% cases were below 5 years of age while >15 years and 5-15 
years age groups represented 33% and 17% failure cases 
respectively. The dhfr double mutant ACNRNI (with mutations 
at 59 and 108 codon) with the highest frequency in the adult age 
group was associated with CQ_ treatment failure. However, in case 
of SP treatment failure, the same haplotype was found associated 
with treatment failure in adult cases (33.3%) as well as under 5 
years of age (16.7%) while haplotype ACIRNL (with mutations at 
51, 59, 108 and 164 codon) was solely associated with treatment 
failure cases in <5 years of age (Figure 6). Overall, the results were 
indicative of the presence of multi drug resistant strains of Pf in the 
study sites. 



Table 2. Prevalence of dhps and dhfr mutations In isolates of Plasmodium falciparum. 





Locus Codons 


Status 


Total (n = 173) 


Location 














Hailakandi (n 


= 103) 


North Lakliimpur (n = 


70) 


dhps 436 


Wild (S) 


130 (75.1) 


70 (68.0) 




60 (85.7) 






Mutant (F/A) 


43 (24.9) 


33 (32.0) 




10 (14.3) 




437 


Wild (A) 


70 (40.5) 


35 (34.0) 




35 (50.0) 






Mutant (G) 


103 (59.5) 


68 (66.0) 




35 (50.0) 




540 


Wild (K) 


106 (61.3) 


48 (46.6) 




58 (82.9) 






Mutant (E) 


67 (38.7) 


55 (53.4) 




12 (17.1) 




581 


Wild (A) 


1 36 (78.6) 


78 (75.7) 




58 (82.9) 






Mutant (G) 


37 (21.4) 


25 (24.3) 




12 (17.1) 




613 


Wild (A) 


134 (77.5) 


75 (72.8) 




59 (84.3) 






Mutant (S/T) 


39 (22.5) 


28 (27.2) 




11 (15.7) 




Locus Codons 


Status 


Total (n = 168) 


Location 














Hailakandi (n 


= 97) 


North Lakhimpur (n = 


71) 


dhfr 16 


Wild (A) 


166 (98.8) 


95 (97.9) 




71 (100.0) 






Mutant (V) 


2 (1.2) 


2(2.1) 




0 




50 


Wild (C) 


168 (100.0) 


97 (100.0) 




71 (100.0) 






Mutant (R) 


0 


0 




0 




51 


Wild (N) 


1 34 (79.8) 


77 (79.4) 




57 (80.3) 






Mutant (1) 


34 (20.2) 


20 (20.6) 




14 (19.7) 




59 


Wild (C) 


5 (3.0) 


2 (2.1) 




3 (4.2) 






Mutant (R) 


163 (97.0) 


95 (97.9) 




68 (95.8) 




108 


Wild (S) 


0 


0 




0 






Mutant (N) 


168 (100.0) 


97 (100.0) 




71 (100.0) 




164 


Wild (1) 


144 (85.7) 


89 (91.8) 




55 (77.5) 






Mutant (L) 


24 (14.3) 


8 (8.2) 




16 (12.5) 




Values within bracket are the percentage of occurrences. 
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dhps dhfr Codons 

Figure 4. Prevalence of dhps and dhfr point mutations in the study. 

doi:1 0.1 371 /journal.pone.01 05562.g004 

Discussion tives, anti-malarial resistance to P. falciparum has been reported 

in vivo for all drugs [30]. Chloroquine resistance to falciparum in 
Control of drug resistant P. falciparum malaria has become a South-east Asia originated from a single focus in Paillang region 
challenging task globally. Except for artemisinin and its deriva- 

Table 3. The frequency distribution of SNPs combination for dhps alleles. 



dhps haplotypes 


Total (n = 173) 


Location 








Hallakandi (103) 


North Lakhlmpur (70) 


SAKAA 


40 (23.1) 


15 (14.6) 


25 (35.7) 


SAKAS/T 


7 (4.0) 


3 (2.9) 


4 (5.7) 


SAKGA 


7 (4.0) 


6 (5.8) 


1 (1.4) 


SAKGS/T 


1 (0.6) 


1 (1.0) 


0 


SAEAA 


6 (3.5) 


4 (3.9) 


2 (2.9) 


SAEAS/T 


2(1.2) 


1 (1.0) 


1 (1.4) 


SGKAA 


1 7 (9.8) 


5 (4.9) 


12 (17.1) 


SGKAS/T 


3 (1.7) 


2 (1.9) 


1 (1.4) 


SGKGA 


13 (7.5) 


6 (5.8) 


7 (10.0) 


SGKGS/T 


9 (5.2) 


6 (5.8) 


3 (4.3) 


SGEAA 


16 (9.2) 


13 (12.6) 


3 (4.3) 


SGEAS/T 


4 (2.3) 


4 (3.9) 


0 


SGEGA 


3 (1.7) 


3 (2.9) 


0 


SGEGS/T 


2 (1.2) 


1 (1.0) 


1 (1.4) 


F/AAKAA 


2 (1.2) 


0 


2 (2.9) 


F/AAEAA 


4 (2.3) 


4 (3.9) 


0 


F/AAEAS/T 


1 (0.6) 


1 (1.0) 


0 


F/AGKAA 


5 (2.9) 


2 (1.9) 


3 (4.3) 


F/AGKAS/T 


2 (1.2) 


2 (1.9) 


0 


F/AGEAA 


20 (11.6) 


16 (15.5) 


4 (5.7) 


F/AGEAS/T 


7 (4.0) 


6 (5.8) 


1 (1.4) 


F/AGEGA 


1 (0.6) 


1 (1.0) 


0 


F/AGEGS/T 


1 (0.6) 


1 (1.0) 


0 



Allelic combinations are in order of S436F/A, A437G, K540E, A581G, and A613S/T, where bold and underlined alleles denotes mutations. Values within bracket are the 

percentage of occurrences. 

doi:l 0.1 371/journal.pone.C105562.t003 
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Table 4. The frequency distribution of SNPs combination of dhfr alleles. 





dhfr haplotypes 


Total (n = ie8) 


Location 








Hallakandl (97) 


North Lakhlmpur (71) 


ACNCNI 


3 (1.8) 


2 (2.1) 


1 (1.4) 


ACNRNI 


108 (64.3) 


67 (69.1) 


41 (57.7) 


VCNRNI 


2 (1.2) 


2 (2.1) 


0 


ACIRNI 


31 (18.5) 


18 (18.6) 


13 (18.3) 


ACNCNL 


2 (1.2) 


0 


2 (2.8) 


ACNRNL 


19 (11.3) 


6 (6.2) 


13 (18.3) 


ACIRNL 


3 (1.8) 


2 (2.1) 


1 (1.4) 



Allelic combinations are in order of A16V, C50R N51I, C59R, S108N and I164L where bold and underlined allele denotes mutations. Values within bracket are the 

percentage of occurrences. 

doi:l 0.1 371 /journal.pone.Ol 05562.t004 



near Thai-Cambodia border during 1950s. Resistance to SP was markers such as pfirt, pfmdr-\, dhfr, dhps & pfalp-6 associated 

detected in 1960s again from Thai-Cambodian border from where with anti-malarial treatment failure provides sufficient evidence to 

Pf strains resistant not only to CQ, but also to SP, mefloquine and policy makers for a renewed thinking on the existing drug policy 

quinine spread to surrounding countries [4]. This spread of anti- [31]. 

malarial resistance in South-east Asia apparentiy followed a Geno typing of circulating Pf strains for C(^ resistance revealed 

pattern until 1978. Initially, it remained limited to tropical forested that almost all isolates were resistant at pfcrt locus, while 68% 

areas corresponding well with the geographical distribution of isolates showed resistance at pfmdr- 1 locus. Simultaneous in vivo 

Anopheles dims complex mosquitoes; thereafter, it moved out of therapeutic study recorded 80% treatment failure to CQ_. Age 

the geographical limits oi An. dims. In India too, the same trend stratified analysis revealed that of all the CQ, failure cases (Jl fii), 

was witnessed where Pf strains resistant to CC^first appeared in 78% (25/32) belonged to <15 years of age including 100% cases 

1973 [5] and to SP in 1987 [6] in An. dims complex mosquito (7/7) belonging to <5 years indicating the synergistic role of host 

dominated areas of north-east India before moving to other parts immunity in the clearance of parasitaemia. Although treatment 

of the country. Therefore, north-east India may be acting as the failure in the older age group (> 1 5 years) was also pretty high 

gateway to resistant falciparum strains in India. Increasing failure (22%), it might partly be explained based on the fact that many 

rates of anti-malarials in north-east India demands a close adult subjects included in the study were not the permanent 

monitoring of the epidemiology and dynamics of the drug resistant inhabitants of the study sites and had migrated from non/less 

Pf Availability of easy and rapid molecular markers would greatly endemic areas for temporary contractual works in search of 

facilitate this process and allows overcoming difficulties in the use livelihood and possibly had low immune status against malaria, 

of traditional methods for characterizing anti-malarial drug Earlier studies in north-east India [2,10] have also reported CQ_ 

resistance. Moreover, the molecular surveillance of drug resistance resistance in more than 80% isolates of Pf indicating non- 

■ HLK(%) ■NLK(%) 




Alleleic combinations 



Figure 5. Frequency of allelic combinations mostly conferring resistance to SP. 

doi:1 0.1 371/journal.pone.01 05562.g005 
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Figure 6. Age and codon wise association of CQ (above) and SP (below) ETF. 

doi:1 0.1 371/journal.pone.01 05562.g006 



effectiveness of CQ, over the region. Althougli /j/mrfr- 1 has been 
thought to affect CQ_ resistance, however, no correlation between 
the prevalence oipfcrt mutation and pfmdr-l point mutation was 
noted in the present study, which was also observed in earlier 
studies further strengthening the doubt over its involvement in CQ, 
resistance [32,33]. 

In case of dhps, though the prevalence of sensitive haplotype 
(SAKAA) was as high as 23%, double (25.4%), triple (22.5%) and 
quadruple (5.8%) mutant haplotypes were also present in good 
numbers along with a full box mutation (F/AGEGS/T) in one 
isolate from HaUakandi district. To the best of our knowledge this 
is the first study recording as many as 23 unique dhps haplotypes 
and a full box mutated haplotype anywhere in India. Earlier 
studies found 8 [34] and 15 [35] dhps haplotypes. Presence of the 
mutated codon 437G (frequency 60%) in most of our dhps 
haplotypes suggests that mutation in this locus started settiing first 
at 437 position in north-east India, similar to the inference made 
by others [26] . 

At dhfr locus, codon 108 was found mutated in all the isolates 
while no mutation was seen in codon 50. Double, triple and 
quadruple mutations were noted resulting in 7 different but unique 
haplotypes of dhfr. Elsewhere in India, presence of double mutants 
(59R/108N) or single mutants (108N) were found to be common 



[22,36]. Although the presence of triple mutants viz. 51I/59R/ 
1 08N or 59R/ 1 08N/ 1 64L were observed from north-east India in 
~10% isolates [22,37], no study has ever reported quadruple 
mutations (51I/59R/108N/164L) from this part of country. We 
found 3 isolates (~2%) harboring quadruple mutants (51I/59R/ 
108N/164L) in the present study (Table 4) which was earlier 
recorded only from Car Nicobar Island, India [37,38] . We found 
the presence of double mutants of dhfr i.e., C59R/S108N or 
N51I/S108N associated with intermediate level of resistance to 
SP, whUe presence of triple mutant (N51I/C59R/S108N) was 
seen in 32% isolates of early treatment failure of SP in therapeutic 
assessment. 

This study revealed high frequency of mutant codons at both 
dhps and dhfr loci that correlated well with the status of 
therapeutic response to SP in Pf malaria cases. Presence of 
haplotypes with mutations at codon position 437 and 540 in the 
dhps gene combined with triple and quadruple mutations in dhfr 
gene, seen in about 10% of isolates comprising of 8 haplotypes 
(Table 4), was predictive of early treatment failure of SP. Overall 
prevalence of mutated codons of dhps and dhfr in north-east India 
detected in our study matches well with the data available in the 
World Wide Antimalarial Resistance Network Molecular Surveyor 
(WWARN: http://www.wwarn.org/surveyor/) (Figure S3). From 
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Table 6. dhps/dhfr haplotypes associated with CQ and SP early treatment failure (ETF). 



dhps/dhfr haplotypes in CQ ETF 




dhps/dhfr haplotypes in SP ETF 




SI. No 


Haplotypes 


Frequency 


SI. No. Haplotypes 


Frequency 


1 


F/AGEAS/T-ACNRNI 


2 


1 SGKGA-ACNRNI 


2 


2 


F/AGEAA-ACIRNI 


1 


2 SGKAA-ACNRNI 


1 


3 


F/AGEAA-ACNRNI 


2 


3 SAKGA-ACIRNI 


1 


4 


F/AGKAA-ACIRNI 




4 SGKGS/T-ACrRNL 


1 


5 


F/AGKAA-ACNRNI 




5 SGKAS/T-ACIRNL 


1 


6 


F/AAEAS/T-ACNRNI 








7 


F/AAEAA-ACNRNL 








8 


SGEGS/T-ACIRNL 








9 


SGEAA-ACIRNI 








10 


SGEAA-ACNCNL 








11 


SGEAA-ACNCNI 








12 


SGKGS/T-ACNRNI 








13 


SGKGA-ACIRNI 








14 


SGKGA-ACNRNI 








15 


SGKAS/T-ACNRNI 








16 


SGKAA-ACNRNL 








1 7 










18 


SAEAA-ACNRNI 








19 


SAKGS/T-ACNRNI 








20 


SAKGA-ACNRNI 








21 


SAKAS/T-VCNRNI 








22 


SAKAS/T-ACIRNI 








23 


SAKAS/T-ACNRNI 








24 


SAKAA-ACIRNI 


4 






25 


SAKAA-ACNRNI 


2 







Allelic combinations are in order of S436F/A, A437G, K540E, A581G, and A613S/T (dfips)-A16V, C50R N51I, C59R, S108N and I164L [dhfr) where bold alleles denotes 
mutations. 

doi:l 0.1 371 /journal.pone.01 05562.t006 



the distribution of the mutations, it can also be inferred that most 
of the mutations in dhps and dhfr loci are fixed in north-east India 
and is now spreading to other parts of India. 

The present study provides important information not only on 
the resistance pattern of CQ, and SP drugs to Pf in Assam state of 
north-east India (2006-2007) but also raises concern on the 
potential effectiveness of the current first line of the drug i.e. 
Artemisinin based combination therapy (ACT) along with 
highlighting the importance of constant molecular surveillance of 
drug resistance. ACT (Artesunate with SP) has been introduced in 
north-east India since 2008 as the first line of drug for treating 
uncomplicated P. falciparum cases. Present study indicates that 
although SP is therapeutically effective in treating uncomplicated 
Pf malaria cases in the study sites, yet fixation and prevalence of 
double, triple and even quadruple mutations at both the loci 
associated with its treatment failure was quite high. 

Although we didn't explore pfatp6, the molecular marker for 
artemisinin resistance, in this study, as ACT was not in use during 
the study period, yet it can be speculated from the findings that 
due to wide spread occurrence of rfA/r and dhps mutant codons in 
north- east India, the artesunate, that is combined with 
sulphadoxine-pyrimethamine in ACT, is probably acting as a 
monotherapy. Since artemisinin is a short acting drug with a rapid 
parasite reduction and a shorter half-life (in hours only), it is 



generally combined with a partner drug (such as SP) having longer 
half life with cUflferent mechanism of action. However, in areas of 
high transmission, like in north-east India, the partner drug may 
promote the proliferation of drug resistant parasites exposed to sub 
therapeutic levels of drug, which may result in decreased efficacy 
of ACT [39]. It is further speculated that the situation in north- 
east India might further worsen by changes in frequencies of 
mutant alleles with the passage of time due to the selective drug 
pressure which raises serious concern over using SP as a partner 
drug with artemisinin for the treatment of Pf cases. Therefore, 
current therapeutic efficacy studies of ACT in north-east India are 
urgently warranted. Although there is a paucity of data on the 
therapeutic efficacy of ACT from this part of India, yet a study on 
ACT involving 4 sites in the Myanmar bordering areas of north- 
east India carried out during 2004—2006 (long before ACT was 
introduced as a first line of drug for treating Pf cases in north-east 
India) recorded 95.5% ACPR of ACT with mean parasite 
clearance time of 55: 17± 14:26 hrs (RMRC unpublished data). 
It was significant that the study revealed 4.5% overall failure (Late 
clinical failure 2.1% and Late parasitological failure 2.4%). This 
might be due to the resistance of circulating Pf strain to SP, since 
most of the treatment failures in ACT is related to the failure of SP 
[40]. In the same study delayed parasite clearance time of 
69:49±6:58 hrs was observed at one of study sites indicating 
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relatively slower effectiveness of ACT in that area. Though 
delayed parasite clearance time is regarded as the most important 
surrogate signal for artemisinin resistance [41] in the absence of 
pharmaco-kinetic study, it is difficult to pin point whether 
increased parasite clearance time observed in that study was due 
to the resistance to artesunate or SP. However, more recently (in 
2013) the Nationsd Vector Borne Disease Control Programme, 
Government of India has changed its drug policy to ACT-AL co- 
formulated tablet of artemether - Lumefantrine combination for 
three days regimen for the treatment of Pf cases of north-east India 
only [42] whereas in the rest of the country a status - quo of using 
ACT (artesunate and SP) for three days as before was maintained. 

A very high level of resistance to CQ,and growing resistance to 
SP, as evident in the present study through therapeutic evaluation 
and molecular characterization of circulating Pf strains, is a serious 
concern for malaria control in north-east India. In view of 
increasing resistance to SP, the partner drug in the ACT 
combination th("rapy used in north-east India, and considering 
this region as a gateway for the drug resistant Pf strains from 
South-east Asia to the South Asia, a close monitoring on the anti- 
malarial drug resistance status is necessary with a constant and 
area wide surveillance for the therapeutic efficacy study along with 
screening of additional molecular markers for artemisinin and 
multidrug resistance (such a.s pfalp6, pfmdrl D1042N and pfmdr\ 
Y1246D etc.,) for overall assessment of the dynamics of 
therapeutic failure of anti-malarials and longitudinal resistance 
monitoring in north-east India to treat Pf malaria cases. Again the 
availability of the recently introduced new ACT-AL combination 
drugs and proper implementation of the new drug policy in NE 
India needs to be monitored from time to time to make the change 
more eflfective and to reduce the emergence of drug resistance as 
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